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Abstract 

We study the effect of quantum ffuctuations on the multiphase point of the 
Heisenberg model with first- and second-neighbor competing interactions and 
strong uniaxial spin anisotropy D. By studying the structure of perturbation 
theory we show that the multiphase degeneracy which exists for 5 = oo (i.e., 
for the ANNNI model) is lifted and that the effect of quantum fluctuations 
is to stabilize a sequence of phases of wavelength 4,6,8,... . This sequence is 
probably an infinite one. We also show that quantum fluctuations can mediate 
an infinite sequence of layering transitions through which an interface can 
unbind from a wall. 

PACS numbers: 75.30.Et, 71.70.Ej, 75.30.Gw 
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I. INTRODUCTION 



There are many naturally occurring examples of uniaxially modulated structures. The 
ferrimagnetic states of the rare earths [ |l| include several phases where the wavevector lies 
along the c-axis and can be of a long period commensurate or incommensurate with the 
underlying lattice. Modulated atomic ordering has been observed in metallic alloys such as 
TiAla and a relationship established between the wavelength of the modulated phases and 
the temperature[ |^. Polytypism describes the phenomenon whereby a compound can have 
modulated structural order of different periods [ ^. A well-known example is SiC where 
the 'ABC stacking sequence of the close-packed layers can correspond to many varied and 
often very long wavelengths. These systems have been usefully modeled in terms of arrays of 
interacting domain walls [ 0. When the wall energy is small, wall-wall interactions become 
important in determining the wall spacing and small changes in the external parameters can 
lead to many different modulated phases becoming stable. 

A model which has proved very useful for understanding this process is the axial next- 
nearest neighbor Ising or ANNNI model which is an Ising system with first- and second- 
neighbor competing interactions along one lattice direction [ |^. At zero temperature the 
ANNNI model has a multiphase point where an infinite number of phases are degenerate 
corresponding to zero domain wall energy. At low temperature entropic fluctuations cause 
domain wall interactions which stabilize a sequence of modulated structures Our aim 

in this paper is to investigate whether quantum fluctuations can play a similar role. That 
quantum fluctuations can remove ground state degeneracies not required by symmetry was 
pointed out by Shender and termed "ground state selection" by Henley [ |10|. 

We find that quantum fluctuations do indeed remove the infinite degeneracy of the multi- 
phase point of the ANNNI model. A sequence of first order transitions is stabilized in a way 
qualitatively similar to the finite temperature behavior but involving a different sequence 
of phases. However, for long-period phases entropic and quantum fluctuations behave in a 
subtly different way. 
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Our analysis focuses on the domain wall interactions and we calculate in turn the wall 
energy, two- wall interactions and three- wall interactions [ This is done by an analysis of 
the structure of perturbation theory around the multiphase point of the ANNNI model: all 
orders of perturbation theory are important. The calculation is described in Sees. 3 and 4 
and corrections pertinent to the long-period phases are treated in Sec. 5. 

To illustrate the essence of the phenomenon, we start, in Sec. 2, by focussing on a simpler 
problem concerning the unbinding of a single interface. In this model, which is effectively a 
one-wall version of the ANNNI problem, spins at opposite sides of the system are fixed to 
be antiparallel. When the magnetic field, /i, is nonzero, the domain wall separating up spins 
from down spins is bound to one of the surfaces. For h = and for an Ising model, there is a 
multiphase degeneracy, because the interface energy is independent of its distance from the 
surface. However, when the Ising model is replaced by a very anisotropic Heisenberg model, 
then, as we show, quantum fluctuations induce a surface-interface repulsion resulting in the 
interface's unbinding through a series of first order layering transitions This calculation 
is similar in spirit, but much simpler than that considered in the rest of the paper. 



II. INTERFACE UNBINDING TRANSITION 

Our first aim is to show how quantum fluctuations can affect the unbinding transition of 
an interface from a surface. Accordingly, we consider the Hamiltonian 

T N-l 1 N-l N TT 

H = E S. ■ S,^i + ;5 E -^^^[31? - S') - -{SI - S^), (1) 

where i are the sites of a one-dimensional lattice of length and Sj is a quantum spin 
of magnitude S at site i. In Eq. (|l|) we introduced factors of S to simplify the classical 
spin (5* oo) limit. Although the results are described for one dimension, they hold for 
any dimension because of the translational invariance of the interface parallel to the surface 
(walls are flat in two or more dimensions for an Ising model at sufficiently low temperature). 
The final term is chosen to impose the boundary conditions such that there is an interface in 
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the system. The interface will be defined as being in position k when it lies between sites k 
and k + 1. We shall restrict ourselves to the limits of zero temperature, H = oo and = oo. 

For D = oo, = aiS where at = ±1 and the Hamiltonian (|I|) reduces to an Ising model 
in a magnetic field whose Hamiltonian Tij is given by 



N-l N-1 

Hi = - J aiGi+i + h'Y (Ti- H[ai - (Jat] . 



(2) 

i=l i=2 

For h > the interface lies at /c = 1; for /i < it unbinds toA; = oo. h = is a multiphase 
point where every interface position has the same energy. For classical spins, S = oo, 
the ground state and hence the multiphase point are maintained as the spin anisotropy is 
decreased from D = oo. 

Our aim here is to study the way in which this degeneracy is lifted by quantum fiuctu- 
ations when D ^ J and S is large but finite. We find that the interface unbinds through 
an infinite sequence of first order transitions as h 0+, as illustrated schematically in 
Fig. 0. The existence of the transitions follows from considering the structure of degener- 
ate perturbation theory around the multiphase point. To start the analysis we write the 
Hamiltonian (|1]) in bosonic form using the Dyson-Maleev transformation [ p!3|JI^ 



S- = (Ji{S - aftti 



25 S.., 



2S 



1 - 



(3) 



where Sa,b is unity if a = 6 and is zero otherwise, af (ai) creates (destroys) a spin excitation 
at site i, and cTj specifies the sign of the ith spin. The resulting Hamiltonian is 



(4) 



where Tij is given in Eq. 



N-l 



^o=E 



i=2 



2D + Jcrj(o-i_i + CTj+i) - h(Ti) 



Q—l + 



(5) 



V4 represents the four operator terms proportional to l/S*^, and V\\ (Vy/) is the interaction 
between spins which are parallel (antiparallel) 



N-l 

V\\ = - JS'^{afai+i + af_^^ai), (6) 

i=2;i^k 

= -JS^^{alal^^ + Ofc+iafc). (7) 

We work to lowest order in 1/5 and therefore neglect terms like V4 which are higher 
order than quadratic in the boson operators. 

To understand the structure of the phase diagram near the multiphase point it is most 
convenient to calculate the energy difference Ai^^ = — E^^i where Ej^ is the energy of 



the system with the interface at position k [ |T6| . In particular, contributions to Ek which 



are independent of k do not affect the location of the interface and need not be considered. 
The energies E^ will be calculated at /i = using standard perturbation techniques [ [T^ 

Ek = k{Om + %)|0)fc - fc(0|(Ml + (Ml + %|0)fc + ■ ■ ■ (8) 

/to — -C/0 

where the vector \0)k corresponds to the configuration with the interface at position k and 
no excitation present and Qo = 1 — |0)fcfe(0|. All the vectors \0)k are eigenstates of Hq with 
the same eigenvalue Eq. However, the perturbative term (y\\ + V[f) conserves J2i Si and thus 
it can never cause a transition between two different ground states. Therefore we may use 
non-degenerate perturbation theory to check whether the excitations can lift the degeneracy 
of the interface states. 

Contributions to the energies Ek arise from spin deviations at the interface created by 
Vy/ which are propagated away from and then back to the interface by V\\ and subsequently 
destroyed by V\j^. However only such processes which are fc-dependent are of interest to 
us. The lowest order term which contributes to AE^ corresponds to an excitation which is 
created at the interface at position k and propagates to the surface and back before being 
destroyed. This graph is illustrated in Fig. |^. (This process contributes to E^, but does 
not occur for Ek-i.) It has a contribution which follows immediately from {2ky^ order 
perturbation theory as 

j2k 1 
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where the terms in J and h in the denominator contribute only to higher order in 1/D. AEk 
is negative corresponding to a repulsive interaction between the interface and the surface 
and hence as — > 0^ the interface unbinds through a series of first order phase transitions 
with boundaries between the phases at 

j2k 

One feature of this calculation which is seen again for the ANNNI model is the fact that 
the interface energy (here AEk) involves the (2/c)th power of the coupling constant, and 
not just the A;th power, as one might imagine for a classsical system [ |T6[. The point is 
that the quantum fiuctuation has to propagate from the interface to the surface AND back. 
As we will see later, this difference leads to a crucial distinction between the way quantum 
fluctuations and classical fluctuations lift the multiphase degeneracy for the ANNNI model. 



III. THE ANNNI MODEL 

A similar formalism is now used to approach a more complicated problem: the effect 
of quantum fluctuations where the multiphase point is a point of inflnite degeneracy for 
bulk rather than interface phases. We take as our example the axial next-nearest-neighbor 
Ising or ANNNI model [ 0. Rather than a single wall interacting with a surface the phase 
structure is now controlled by an inflnite number of interacting walls and we shall follow 
Fisher and Szpilka [ ^ in analyzing the phase structure in terms of the interactions between 
the walls. A brief account of this work has been published elsewhere [ ITS . 

The Hamiltonian we consider is 

^ = X! ■ Sij' - Sjj ■ Sj+ij + t;^ X! ^iJ ' ^*+2j - ^{[^ijT ~ "S"^)' (11) 

^{n ) *J *J 

where i labels the planes of a cubic lattice perpendicular to the 2;-direction and j the position 
within the plane. Also {jj') indicates a sum over pairs of nearest neighbors in the same plane 
and Si J is a quantum spin of magnitude S at site {i, j). For D = oc, only the states S'f = aiS, 
where cxj = ±1 are relevant and Ti reduces to the ANNNI model [ ^. 
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(12) 



The ground state of the ANNNI model is ferromagnetic for k = J2/J1 < 1/2 and an antiphase 
structure with layers ordering in the sequence {cTj} = {...1,1, —1, —1, 1, 1, —1, — 1 . . .} for 
K > 1/2. K = 1/2 is a multiphase point[||, |^, where the ground state is infinitely degenerate 
with all possible configurations of ferromagnetic and antiphase orderings having equal energy. 
For classical spins, 5* = 00, the ground state (and therefore a multiphase point) is maintained 
as D is reduced from infinity. 

To describe how the degeneracy is broken at the multiphase point when S is large, but 
not infinite, we define a notation similar to that of Fisher and Selke[ |^ using {ni,n2, ■ ■ ■ fim) 
to denote a state consisting of domains of parallel spins with alternate orientation whose 
widths repeat periodically the sequence {rii, n2, ... rim}. 

As in the previous section we use the Dyson-Maleev [ [T^,[I^ transformation to recast the 
Hamiltonian ([TT]) into bosonic form (working to lowest order in 1/S) with the result 



ni{a,}) = Eo + no + Vii + V^ 



(13) 



where Eq = Ha, 



Ho = E 



2D + Ji(Jij((Ji_ij + CTj+ij) - J2<^i,j{ai-2,j + (Tj+2,, 



c-l + 



(14) 



with D = D + 2Jq and V\\ (V^) is the interactions between spins which are parallel (antipar- 
allel) 



-JiX{i,i + l;j){afjai+ij + af^-^^aij) + J2X{i,i + 2; j)(a+-ai+2,i + «i+2,j«i,i) 

(15) 



-JiY{i,i + 1; j)(ajaj^ij + aj+ijajj) + J2Y{i,i + 2; j)(a+.a+ 2^ + ai+2,jC 



(16) 



where X{i,i';j) \Y{i,i';j)] is unity if spins and are parallel [antiparallel] and is 

zero otherwise. We do not consider quantum fluctuations within a plane, since the phase 
diagram is determined by the interplanar quantum couplings. Moreover we shall work to 
leading order in 1/5* , in which case four-operator terms can be neglected [ |T5|. Also we will 
continue to use non-degenerate perturbation theory, since the perturbative term (y\\ + Vy/) 
cannot connect states in which the wall is at different locations, since such states have 
different values of J2i S^. 

The structure of the phase diagram will be constructed by considering in turn E^,., the 
energy of an isolated wall; V2(n), the interaction energy of two walls separated by n sites; and 
generally Vfc(ni, n2, . . . n^.i), the interaction energy of k walls with successive separations 
rii, n2, ... Uk^i [ If. In terms of these quantities one may write the total energy of the system 
when there are walls at positions rrii as 

E = Eo + n^E^ + ^2("^i+i -m,) + Yj ^3("^i+2 - ^^+1,771^+1 - rrii) 

i i 

+ V4(mj+3 - mi+2, mi+2 - m+i - nii) + . . . , (17) 

i 

where Eq is the energy with no walls present. The scheme of Ref. [ for calculating the 
general wall potentials Vk is illustrated in Fig. |^. Let all spins to the left of the first wall 
have (Jj = 0" and those to the right of the last wall have ai = rj for k even and (Xj = —rj for 
k odd. The energy of such a configuration is denoted Ek{cr,rj). If cr = — 1 {rj = —1) the left 
(right) wall is absent. Thus the energy ascribed to the existence of k walls is given by[ |19[ 

Vfc(ni,n2, . . .nfe_i) = Yl <yr]Ekicr,r]) . (18) 

a,ri=±l 

Contributions to E/^ which are independent of a or do not influence V^. Ek{a,T]) is cal- 
culated by developing the energy in powers of the perturbations Vy/ which allows creation 
(and annihilation) of a pair of excitations straddling a wall and V\\ which allows the excita- 
tions to hop within domains. We consider contributions to the wall energy and to two- and 
three-wall interactions in turn. 
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A. Wall energy 



Contributions to the wall energy to second order in perturbation theory arise from ex- 
citations which are created at a wall and then immediately destroyed as shown in Fig. ^. 
These effectively count the number of walls and therefore lead to a renormalization of the 
wall energy of 

ADS \D^SJ ^ ' 

but since we work to leading order in S*"^, the 5*^^ correction to will not affect the results 
for Vk. 



B. Pair interactions 

The lowest order contributions to V2(n) are obtained by creating an excitation at, say, 
the left wall using and then using V\\ for it to hop to the right wall and back. Because 
we assume the existence of the left wall, this contribution implicitly includes a factor d^^i. 
Now we look for the lowest-order (in J / D) contribution which also has a dependence on rj. 
In analogy with the unbinding problem, we might consider processes in which the excitation 
hops beyond the wall. Since such a term can not occur when the wall is actually present, 
it will carry a factor For n odd, we illustrate this process in Fig. ^, and see that it 

gives a contribution to V2{n) of order Jg^^/D". As we shall see, there is actually a slightly 
different process which comes in at one order lower in J/D. To sense the presence of the 
right-hand wall, note that Eij in Eq. (|1^ will depend on r] if the i is within two sites of the 
wall. Therefore it is only necessary to hop to within two sites of the right wall, as shown 
in Fig. ^ for an energy denominator {Hq — Eq) in the series expansion (|^) to depend on rj. 
This process is of lower order in J/D because it takes two interactions to hop back and forth 
but only one to sense the potential via an energy denominator. Accordingly, in contrast to 
the interface unbinding considered in Sec. 2, it is necessary to retain the terms in the J's 
in the energy denominators to obtain the leading order contribution to V2{n). We consider 
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separately n odd and n even. 
n odd 

To lowest order the processes which contribute are those shown in Fig. |^. For a domain of 
n spins with cTj = — 1, (n — 1)**^ order perturbation theory gives 

E2{c7,r^) = 25„^iJ^-^S-\-lY~^{Ab + 2Ji}-2{4D + 2Ji - 2J2}-("-5) 

x{4L) + 2Ji- J2(l-r/)}-^ (20) 

In writing this result we dropped all lower-order terms because they do not depend on both 
a and rj. Here and below, the dependence on a is contained in the factor because we 
assume the existence of the left-hand wall. The energy denominators are constructed as 
follows. The left-hand excitation has energy 2D since it is next to a wall. The right-hand 
excitation has the energy according to its position as illustrated in Fig. ^ The prefactor of 
2 arises because the initial excitation can be near either wall and the overall factor (— l)""^ 
arises from the (—1) associated with each energy denominator. Adding the contributions 
from (|20|) appropriately weighted as in ([T8|) gives 



{AD + 2JiY{AD + 2 Ji - 2 Js}"-^ L 4/) + 2 Ji AD + 2 Ji - 2 J; 



= 4J2"^-V(4Z})"-i + C(l/Z}"), nodd. (22) 

Note that there is no term (9(1//)"^^). This is because to this order the energy denominators 
are independent of the J's. Hence to this order Ek{a, rf) is independent of 77 and the sum in 
Eq. (0) is zero. Similarly terms from n*'* order perturbation theory (in which one J2 hop is 
replaced by two Ji hops) do not contribute C(l/I)"^^). 
n even 

For even n several diagrams contribute to leading order, i. e., at nth order perturbation 
theory. These are shown in Fig. |^. As an example we give the contributions to the energy 
from the diagram (b)(iii). Again we drop all terms which do not depend on both a and rj. 
Thus 
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2 

(4D + 2Ji - 2J2)-^"-^\4D + 2Ji - J2(l - T])]-' , (23) 

where the superscript iii indicates a contribution from diagram iii of Fig. ^ the prefactor 2 
comes from including the contribution of the mirror image diagram, the prefactor (—1)"^^ 
is the sign of nth order perturbation theory, the factor {n — 2)/2 is the number of places the 
single (Ji) hop can be put, and indicates that this contribution assumes the existence 
of the left-hand wall. To leading order in D, the 77-dependence is contained in 

Et\^v) = {-ir-\n-2)S^,,j',jr'S-\4Dp^-'\4D + 7^de2/dr^)-' 

~ (-l)"(n - 2)r]6,^^J^rr'S-\4Dy''{de2/dr]) . (24) 

Using de2/drj = J 2, we get 

V^f = 2{n - 2)J^J^-^S-\AD)-^ . (25) 

We treat the other diagrams of Fig. |^ similarly. Dropping terms which do not depend on 
both a and 77 and working to lowest order in {D)~^, we get 



2Jl{de2/d7]) + ]^Jl{n - 2f{de2/d7]) 



E2{a,v) = r]6,,,Jr'S-\4D)-^ 

+ {n - 2)Jl{de2/dri) + (n - 2)Jl{de2/ drj) + 2J^{dei/dri) + 2J^{de2/dri) 



(26) 



where the contributions are from each diagram of Fig |^, written in the order in which they 
appear in the figure. Thus for n even we have 



V2{n) = S-^ADyX 

jn—l 
Jo 



n jn—2 
2 



AJfJ2 + JtJ2{n - 2Y + 4(n - 2) + A4{J2 - Ji) + 
n^Jl - U1J2 + 8J^), n even, (27) 



{ADys 

where we used de2 /drj = J2 and dei / drj = J2 — Ji- 

Fisher and Szpilka [ Q have shown that the phase sequences can be determined 
graphically by constructing the lower convex envelope of V2{n) versus n. The points 
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[n, V2(n)] which he on the envelope correspond to stable phases. The pair interactions 
defined by the expressions (^) and ( pTf ) correspond already to a convex function for 
n « (i)/J)^/^. Hence, in this regime, we expect within the two-wall approximation 
a sequence of phases (2), (3), (4), . . . as shown schematically in Fig. ||. The widths of 
the phases (n) can be estimated using the fact that each phase is stable over an interval 
AEu, = n[V2{n~l) — 2V2{n) — V2{n + 1)] [Q. Therefore, using (\[% we can say that the width 
A( J2/ Ji) occupied by the phase (n) in Fig. | is C(( Js/L*)""^) for n odd and 0(( Js//^)""^) 
for n even. This sequence of layering through unitary steps (n) — > (n + l) will not be obeyed 
for large n, i. e., for n ~ (D/jy^'^, because then V2{n) will suffer from strong even-odd oscil- 
lations. Moreover, for large n, the entropy of more complicated perturbations may dominate 
the physics. A discussion of this is given in Sec. 5. Here we go on to consider the effect of 
3-wall interactions which can split the phase boundaries (n) : (n + l) where there is still a 
multiphase degeneracy of all states comprising domains of length n and n + l. 



IV. THREE- WALL INTERACTIONS 

Three-wall interactions are needed to analyze the stability of the (n) : (n + l) phase 
boundary to mixed phases of (n) and (n + l). The condition that the boundary be stable is 

[i 

F{n, n + l) = V3{n, n) - 2V3{n, n + 1) + ¥3(71 + 1, n + 1) < 0. (28) 

Consider first the calculation of F{2n — l,2n). The diagrams which contribute in leading 
order to V3(2n — 1, 2n — 1) and V3(2n, 2n — 1) are shown in Figs, ^a and P?, respectively. 
To leading order in 1/D, V3{n + l,n + 1) does not contribute to F{n,n + 1). Figure |^ 
aims to emphasize the positions of the initial excitation and the closest approaches to the 
neighboring domain walls. One must also consider the position of the first neighbor hops in 
B and C and the sequence of the hops when calculating the contribution of the diagrams. 
An explicit calculation of the contributions of the relevant diagrams would be extremely 
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tedious. However what concerns us here is the sign of F{2n — 1, 2n). If Ni is the contribution 
to F of diagrams of type i in Fig. |^ , 

F{2n - 1, 2n) = 2Na + 2Nb + 2Nc - 2Nq, (29) 

where the factors of 2 multiplying A^a, and A^^c account for the mirror image diagrams 
and that multiplying occurs because of the 2 in Eq. (P^). 

We shall now show that F{2n — 1, 2n) < 0. Consider a diagram in which the hops occur 
in the same order in A, B, C and D and the Ji hops in B and C are, say, nearest the outer 
walls. The matrix elements of all types of diagram carry a negative common factor (the 
sign arising because we are considering even-order perturbation theory) and their ratios are 
rriA/Tn-o = 1 and mxlm^ = mp^/mc = J^/Jf- 

We must also expand the difference in the energy denominators in a way analogous to 
the step between equation ( pT]) and (|2^), but here to second order in J/D. Using ([T8|), the 



contribution of each diagram to the appropriate A^, may be written 



cr,ri 



h H 77^-. ^ 



L(4D)4"-55 (4D) {ADf 



= + Oi^^) (30) 

where the coefficients / depend only on Ji and J2. When the sum is taken only the term 
/lo multiplying err/ survives. For diagrams of type A, /lo is J2(<^2 ~ -^i), while for B, C and 
D it is J|. Therefore these diagrams give a contribution to F proportional to 

- j|Ji(2Ji - J2) < . (31) 

The contributions to F of the other diagrams in B and C (which correspond to a different 
position of the first neighbor hop) is proportional to — J^J|. Hence F{2n — l,2n) < and 
the {2n — 1) : {2n) boundaries are stable. 

A similar argument holds for F{2n, 2n + 1) for n > 1. The relevant diagrams are shown 
in Fig. im They contribute 
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F(2n, 2n + l) = 2Na + 2Nb + 2Nc + 2Nb + 2Ne - 2Nf. (32) 

Using the same argument as above 

A^a + A^b-A^foc-J2Vi(Ji- J2) <0. (33) 

A^'c, Nd, Ne, and the other orderings of A^b are negative and hence F{2n, 2n+l) < 0. Thus 
the phase boundaries {2n) : {2n + 1) are first order for n > 1. 

For the (2) : (3) boundary different diagrams contribute to F(2,3). Indeed the sec- 
ond order expansion of the energy denominators [as in Eq. (pO|)] gives a zero contribution. 
Accordingly, the calculation of F{2, 3) requires going to higher order in {J2/D). This calcu- 
lation is carried out in detail in Appendix A and shows that the (2) : (3) boundary is also 
stable. 



V. LARGE n ANALYSIS 

For small n, we have seen that the leading contribution to V2(n) is of order D{J2/DY/S, 
where the value of x corresponds to the minimum number of steps needed to go from near 
one wall to near the other one and back: x = 2[n/2] + 1, where [x] is the integer part of x. 
As n increases, the contributions from longer paths, although individually less important, 
can become dominant because of their greater entropy. To allow for this possibility we now 
carry out perturbation theory in terms of the exact eigenstates for one excitation in each 
block of parallel spins. In this formulation, the unperturbed Hamiltonian is the sum of 
the Hamiltonians of each domain of parallel spins when all interactions with neighboring 
domains are removed. Thus from equations ([l^) and (^) the unperturbed Hamiltonian for 
a block of parallel spins from sites I to J inclusive can be written 

ni'''^ = Y^JrjS-\at - a+)(a, - a,) + ^2D5-ia+a, , (34) 

id i 

where Jij = Ji5j,i+i — J2'^j,j+2 and the prime on the summation indicates that the sum is 
restricted so that both indices are actually in the block. The matrix representation of Ti^''^^ 
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is given explicitly in Appendix B. It follows from and ([T6|) that the perturbation 
associated with a wall between sites s and s + 1 is 

where 



^S-'J2wt^ata,, (36) 



k 



and = 

In this formulation it is not natural to calculate V2{n) directly. Instead one calculates 
the total energy of given configurations from which ¥2(11) is easily deduced. We start by 
calculating the total energy of a configuration with a single wall between sites and 1. This 
gives the wall energy as 

^(2) ^ (o|y(°)^X(°)|0) = S-'j:xS^xS\0\a,a,^atat\0) (38) 
^ ijki ^ 

where S = Eq — Hq, with Eq the ground state energy, defined to be zero in this context. 

Here we have introduced the notation that the subscript on E specifies the number of walls, 

the superscript the order in perturbation theory, and the arguments (if any) the separations 

between walls. 

To evaluate ( p8[ ) and similar expressions we now introduce the exact eigenstates for a 
single excitation on either side of the wall when interactions across the wall are ignored. 
For a block of parallel spins occupying sites / through J, inclusive, these single-particle 
eigenstates satisfy 

T.{ni!''^)/t'\j) = S-'el^''U'I''\^. (39) 
j 

Later we write e*-^''^^ 
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To evaluate Eq. (|38|) in terms of the exact eigenstates notice that afa^ connects the 
ground state to a state in which the semi-infinite chain to the right of the wall is in an 
excited state which we label (3 and the semi-infinite chain to the left of the wall is in an 
excited state a. Thus we have 

1 Z-^ ij kl (oo) (oo) ■ \ ' 

ijkl a(3 -|- €p 



This process is illustrated in Fig. 11 



We now construct the energy of a system with only two walls, one between sites and 1, 

the other between sites n and n+1. The contribution to the total energy of this configuration 

(1'] 

from second-order perturbation theory, denoted E2 {n) comes from an expression similar to 
Eq. (|38D but which here involves one semi-infinite chain and one block of length n. 

Here and below we include a factor of 2 because the process could be initiated at either 

(2) (2) 

of the two walls. Note that as n ^ 00, E2 2E{ , as expected. We will also need the 
contribution to the energy of this configuration from third-order perturbation theory. The 



only process at this order is shown in Fig. |I2| and it gives a contribution 



(n) =2(0|r(°)|^W^(«)|ix(°)|0) = 25-3 ^ X^Jw^^ X^f {0\a,a,^^ata,^at4\0) . 

ijklm 

(42) 

In Eq. we see that X^^^ creates one excited eigenstate in the semi-infinite chain to the 
left of the walls and also an excited eigenstate in the down-spin block of length n. That 
type of reasoning allows us to rewrite Eq. (^) as 

(n) = 25- E Y.X^fwt'xi'j'^''^'^''^^^^^^^^^ 



^ U k ^•■Im / (00) / (00) 

ijklm a/37 V /3 / I ^'^ ~r C 



7 , 



(43) 



Then, up to third-order perturbation contributions, the wall potential V2{n) we wish to 
obtain is given by 
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V2 {n) = Ef'^ {n) - 2E?^ + E^? {n) . (44) 

To interpret these expressions it is convenient to express them in terms of the Green's 
function, defined by 

Gr(^)^E^^»^. (45) 

a E -\- da 



Thus 



ijkl a 



E?\n) = -25-iEE^.?^4?Vi-~'°H^)0L-°°'°n^)Gj.l-"^(eL~)), (47) 

ijkl a 



and 



ijklm O! 

To obtain V2{n) we will have to determine SG = G^^'"'^ ~G^^'°°\ To evaluate this quantity 
we need to identify the perturbation which, when added to the unperturbed Hamiltonian 
describing two independent blocks of spins, {l,n) and (n + l,oo), gives the unperturbed 
Hamiltonian 7^0^'°°-*. This perturbation V can be written as 

V = + , (49) 

where Z^"-^ describes hopping across the wall which is needed to make the semi-infinite chain 
from a finite block of parallel spins. 

We now use some results of standard perturbation theory for a Green's function, as given. 



for instance, in Ref. 0. For this expansion we work to lowest order in the wall perturbation, 

V of Eq. (^9|). We choose TYq to be the Hamiltonian for a block of n parallel spins and treat 

V perturbatively. In first-order perturbation theory for V, it is not necessary to keep Z*^") 
(and consequently Ti'^'^^''^^) because it moves an excitation to the right of the right wall 
which cannot be hopped back to the (1, n) block without going to higher-order perturbation 
theory. So correct to first order in perturbation theory we have 
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G 



(l,oo) 



-1 



E 



n) 



-1 



in) 



ik 



kj 



where W^^'^ is defined in Eq. (|36|) . Thus using Eqs. 
result 

= AS-' E E^rw^i"^^^S0L-°°'°n^)0L-°°'°nOGS'"He(-))4^^^^ . (si) 

ijklm Oi 

Evaluating this when Ji = 2J2 we obtain [writing G for G^'''^\e'^) and for 0^^°°''^''] 

V^2(n) = ^ E|<^-(0)' (^2,n-l - 2Gi,„_i^ - (G2,n ' 2Gi,„ 



(50) 

8|) and f^Ol) we have the 



+1 



(52) 



4J| 



E< 



n 2 



^'a(O) G2,n-l-2Gi, 



n-1 



-1)G 



l,n-l 



,(0) G'2,„-2G'i,„ ^^^(-l)^,^ 



(53) 



We will evaluate this with successively increasingly accurate approximations for large n. 
For small n it is certainly correct to replace e^°°^ by 2D' = 2D + 2Ji — 2J2 = 2D + 2J2, 
since corrections will be proportional to J/D' with a bounded coefficient. For the moment 
we continue to use this approximation even for large n. With this approximation, the sum 
over a in Eq. (|5T| ) yields 



(54) 



We refer to this as the nonpropagation approximation, since it amounts to setting the off- 
diagonal elements of ?io~°°'°'' to zero, forcing i and I to coincide. 

Within this approximation and writing writing G for G*^^'")(2Z}'), we find that 
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V2{n) = Va + Vb (55) 

where 

4 73 / \ 2 

VA=^[G2,n-l-2Gi,n-l) , (56) 

= ^ (4Gi,„G'i,„_i - bGl)i . (57) 

In Appendix B we give an essentially exact evaluation of the required G's, apart from an 
overall scale factor which we only obtain approximately. This evaluation leads to the result 

V2{n) = ^ ( sin'lnS + 4S'] - (3/8)( J^/I?')'! , n even ; 



^ (^cos2[n5 + - (3/8)( J^/D')') , ^ odd ; 



(58) 



where, to leading order in J2/{4:D'), = = J2/(4D')- Here Va gives rise to the term 
involving the square of the trigonometric function and, when Va = 0, 



\D 



For small n these expressions reduce to our previous results ( P^ and ( P7D at leading order 
in J2/D', in which case Vb is irrelevant. 

We now discuss the interpretation of these results. For the moment let us ignore com- 
pletely the term Vb- When V2{n) is nonmonotonic, as we found here, an elegant graphical 
construction which yields the phase diagram was suggested by Szpilka and Fisher [ Q . This 
proceeds by drawing the lower convex envelope of the points V2{n) versus n. Points on the 
convex envelope are the allowed stable phases (assuming no further bifurcation due to V3). 
For this construction it is important to distinguish the case when V2{n) becomes negative. 
If this occurs, then there will be a first-order transition from rio to = cxd where riQ is the 
value of n for which V2{n) attains its most negative value. On the other hand, if V2{n) is 
positive for all n, then one has an infinite devil's staircase, with no bound on the allowed 
values of n. Accordingly, it is obviously important to ascertain whether or not V2{n) is 
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positive definite. Eq. (|58D suggests that V2{n) can become negative when {nS + 46^) / {27i) is 
sufficiently close to an integer. However the approximations inherent in its derivation may 
alter this conclusion. 



A. Effect of Allowing Propagation of the Left Excitation 

To determine whether or not an unending devil's staircase actually exists in the phase 
diagram, it is necessary to assess the validity of the nonpropagation approximation. We now 
avoid the approximate treatment of Eq. (^Tp in which we replaced e^°°^ by 2D'. We write 
Eq. dlD as 



ijklm c 

where 



a 



I ,Hir,(oo) I (n)i 
t« ' J F^i ~r 67 J 



a [2D' + ef^] [2D' + e^")] V 2D' + ej^^ 2D' + e?^ 



2dD'\2D' + eP 2D' + e^^\ 



= Yo + 6Y . (61) 

Keeping only the term Yq leads to the nonpropagation approximation, and thence to Eq. 
(0) and the results of Eq. (||). 

We now analyze the effect of 6Y. For that purpose we use the fact that the eigenfunctions 
satisfy Eq. (^). For sites i near the wall, (i.e., i = and i = —1) Eq. (|H^) yields [omitting 
the cumbersome superscripts {—00, 0)] 

(e„ - 2D')0„(O) = (J2 - Ji)MO) - Mai-1) + J2</>a(-2) (62) 

(e„ - 2D')0„(-1) = J20«(-l) - Ji0«(O) - Ji0«(-2) + J20„(-3) . (63) 

Using these equations and also Eq. (0), we get 
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5Y 



X 



'I d ( 1 1 



(64) 



5y leads to contributions with i ^ I shown in Fig. O, and also with i = I which are not 
shown but are similar to those of Fig. |^. If G denotes G(^'")(2D'), then we have, from Eq. 
a correction to the two-wall interaction of 



d 



ijklm 



k,m 



X 
2 d 



SdD 



— 4Ji j|(^J2G'2,n-lG'i,„_l + (J2 — Jl)Gi^nG2,fij ~^ '^J^.^J^GX^^i + (J2 — Jl)G\.^j 

+2(J2 - Ji)Jl{j2Gl^_, + (J2 - Ji)G|„^ + 2(J2 - Ji)Jl[j2Gl^_, + (J2 - Ji)Gl)i 



-A{J2- Jl)JlJ2[j2Gl ,n-lG2,n-l + {J2 ~ Jl)Gi^nG2,-. 

We simplify this by setting Ji = 2J2. Then 

2J| 9 



(65) 



5V2{n) 



(66) 



5 9Z}\ 

The dominant contribution comes from the first term. To evaluate this expression, it is 
necessary to develop an expression for G2,n-i- Using Eq. ( [B4a| ) of Appendix B, we write 



A„(-l)"+lG'2,„-l = C2d„_3 = jrl(C/J2) V-'Qn-3 = Jr^y'^^^Qn- 



(67) 

where y = ^ AD' / J2 and we set G = AD'. The calculations for n odd and n even are similar. 
Here we do them only for n even, in which case 



j2n-2 



(68) 



To take the derivative note that Gg^-i ~ Z^' " ^sin^^nS — 26) and 5 ^ D' ^1"^. Thus we 
have 
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= ^y""^ sin^(n5 -26) -2 sin(.5 - 25) cos(n5 - 25)^^^ 



(69) 



Then we obtain, for 1 and n^/vr an integer, 



dD' " A2 ' ^^^^ 



so that 



^''^w = w"^^ = w"^ (!) pi) 

Since ?t,(5 > 1 when Va = 0, this term is larger than Vb by at least {D'y/J2) ~ {^lD' j J'lfl'^ . 
Since this term is positive, we see that allowing for the left-hand excitation to propagate 
leads to a correction to Va which is much more important than Vg. In other words, this 
more accurate evaluation gives V2(n) = Va + 5V2(n). This result relies on the validity of the 
expansion in Eq. (|6l|) , the precise condition for which is not obvious. However, when n gets 
sufficiently large, this expansion breaks down and the considerations in the next subsection 
become necessary. 



B. Large n Limit with Propagation 

The expansion that we have used in Eq. (pT]) implicitly assumes that the Green's function 
has a weak dependence on energy. That is true as long as n is small enough. But when 
n becomes arbitrarily large, then there must exist a regime in which the right-hand side of 
Eq. (|53D is dominated by the largest term in the sum over a. If it were correct to keep only 
a single value of a, then it would be possible to fix D\ so that the first square bracket in 
Eq. (|53D would vanish and V^in) would be negative. This reasoning is not correct, however 
as the analysis in Appendix C shows. Even for large n the sum over a has a width in a 
of order which prevents (G2,n-i — 2G'2,n)^ from being fixed to be precisely zero. We 

therefore conclude that for the one- dimensional system of walls, V2{n) does remain positive 
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for n ^ oo. It seems unlikely that in a crossover between the regimes we have considered 
¥2(11) would become negative. So we conclude that for the one- dimensional problem V2{n) 
remains positive and there is no cutoff in the Devil's staircase for the phase diagram. 

C. Large n Limit For Three-Dimensional Systems 

In the discussion up to now, we have treated the three-dimensional system as if it were 
a one-dimensional system in which planar walls separate up-spin segments from down-spin 
segments. Here we give a brief argument which suggests that these one- dimensional results 
continue to hold for the three-dimensional system. One way to phrase the argument is to 
note that when D' is large compared to the J's, we are far from criticality. The correlation 
length (of order | ln(J/D')|^^) is very short. Thus, entropic effects of longer paths are 
strongly cut off by the correlation length. Here we indicate the nature of a formal argument 
of this type. 

To analyze the three-dimensional case, we consider only the dominant term, illustrated 
in Fig. |n[ It gives rise to the contribution 

4J| 

(n) = — H <^".9±(05 0)<^"-g±(05'^-L + "5-L)G'2,n-i('^±;ea,g^)G'2,„-i(s_L;ea,g^) , (72) 

where we omit the superscripts. The subscripts on (p are the quantum number, a, associated 
with the coordinate perpendicular to the wall and the wavevector associated with the 
transverse coordinates. The arguments of (p are the coordinate perpendicular to the wall and 
the vector displacement in the plane of the wall. The arguments of G are the displacement 
in the plane of the wall and the energy. Considering only the dependence of and e on 
wavevector we obtain 

(5V2(n) ~ exp[iiq_L ■ (r_L + s_L)]G'2,„-i(r_L;eg^)G'2,„_i(s_L;eg^) . (73) 

r±,s±,q± 

In terms of Fourier transformed variables for coordinates in the plane of the wall (indicated 

by overbars), we have 
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SV2{n) - Y^[G2,n-i{q±;eqjf 



(74) 



But this is again the type of expression analyzed in Appendix C. So we conclude that for 
the three-dimensional system V2(n) is also always positive. 



VI. DISCUSSION 

The aim in this paper has been to demonstrate how quantum fluctuations can lead to 
interactions between domain walls and hence stabilize long-period phases in the vicinity of 
a multiphase point where the intrinsic wall energy is small. 

We first considered a Heisenberg model with strong uniaxial spin anisotropy D and an 
interface pinned to a surface by a bulk magnetic field h. A perturbation expansion in D"^ 
was used to show that the wall-interface interaction is repulsive and hence that the interface 
unbinds from the surface through an infinite number of layering transitions as h passes 
through 0. 

The bulk of the paper was devoted to describing the behavior of the Heisenberg model 
with first- and second-neighbor competing interactions and uniaxial anisotropy D near the 
ANNNI model limit D — oo. This model has a multiphase point for sufficiently large D 
which is split by quantum fluctuations to give a sequence of long-period commensurate 
phases (2), (3), (4), ... (n) ... . 

The phase sequence could be established for n not too large by a calculation of two- wall 
and three-wall interactions using perturbation theory with as a small parameter. A 
discussion of correction terms important for large n was given from which we concluded 
that, unhke for the ANNNI model, the sequence of phases is infinite. The reason this model 
is different in this regard from the ANNNI model is an inherently quantum one: for one wall 
to indirectly interact with another an excitation has to propagate from one wall to the other 
AND return. Thus the interaction in the quantum case is proportional to the square of a 
oscillatory Green's function whereas in the ANNNI model the analogous function appears 
linearly. As a consequence of this oscillation, the phases come in the sequence n — > n -|- 1 or 
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n n + 2, depending on the value of n. In the latter case, we did not explicitly investigate 
the stabihty of the phase diagram, but a cursory analysis leads us to believe that the function 
F{n, n + 2) analogous that in Eq. ( pH]) is negative. 

Similar behavior is observed in both the interface model [ [1^ and the ANNNI model 
[ ^, for finite temperatures. Here thermal fluctuations replace quantum fluctuations in 
mediating the domain wall interactions. Although long-period phases are stabilized in the 
ANNNI case, the qualitative form of the phase sequence is very different to that discussed 
in this paper. The stable phases are (2'^3), k = 0, 1,2, ... /cmax with /cmax — > oo as T — 0. 
Mixed phases (2^32^-i3) also appear [ ^ . 

A third mechanism that can split the degeneracy at both a bulk [ ^ and an interface 
2l| ] multiphase point is the softening of the spins themselves; a non-infinite spin anisotropy. 



This does not occur for the ANNNI model where there is a finite energy barrier for the spins 
to move from their positions at D = oo. However, for a similar model with 6- fold anisotropy 
an infinite number of phases become stable near the multiphase point as D is reduced from 
infinity. [ED] 
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APPENDIX A: STABILITY OF THE (2): (3) BOUNDARY 

In general F{2n,2n + 1) is 0{1/D'^^~^). However, for n = 1 the term 0{1/D^) is 
accidentally zero and terms 0{1/D'^) must be retained. This leads to a lengthy calculation. 
We now calculate F{2, 3) explicitly, considering each order of perturbation theory in turn. 

Second-order perturbation theory 

Contributions arise from diagrams spanning a wall which are created and then immedi- 
ately destroyed (as in the example of Fig. 1). The contribution to ¥3(2, 2) comes from both 
first and second neighbor excitations. That to ¥3(2, 3) is just from second neighbors because 
the energy denominator of the first neighbor excitation does not depend on both a and 77. 
For the same reason there is no contribution at all to V3(3, 3). 

Using a subscript 2 to indicate that we are considering only the terms arising from second 
order perturbation theory one obtains 0{1/D^) 

Vi(2, 2) = J^l-Jf + - J2)] + ^^l-J^J^ + ^J^], (Al) 

V,(2.3) = -^ + l?^a±ill, (A2) 

V2(3,3) = 0. (A3) 

Third-order perturbation theory 

Contributions to V3(2, 2) in third order perturbation theory arise from diagrams like that 



shown in Fig. |I5|. Recalling that the spins on either side of the wall can hop, and that the 
initial excitation can be between second neighbors, with a subsequent hop to first neighbors 
gives 

V3(2.2) = l?Me^^. (A4) 
Similar diagrams contribute to V3(2, 3). The hop must lie within the domain of 3 spins 

^3(2,3) = ^^^. (A5) 

There is no contribution V3(3,3). 
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Fourth-order perturbation theory 

We first consider processes which are proposrtional to ^fV"^- As we discussed in the 
text, to lowest order in J2/D we do not need to consider processes which hop beyond the 
wall. However, since the calculation of F(2, 3) requires a calculation of ¥-^{2, 2) and ¥^{2, 3) 
including the first higher-order corrections, we need to keep such processes. 

We now evaluate contributions from such processes, which we show in Fig. First of 
all, since these processes only exist in the absence of the right-hand wall, they all carry a 
factor 5^,-1. Secondly their overall sign is negative for even-order (fourth-order) perturbation 
theory. Also, the contributions to V^{2, 2) carries a factor of 2 to account for the mirror image 
diagrams. 

Thus from the first diagram of Fig. |16|, we get Ek{a, rj) of Eq. (p!8D as 



Esia,r]) = -25r,,-iJlJl{SrS2S^)-^ , (A6) 

where Si is an energy denominator. We have that Si ~ AD + 0(J). In particular, we need 
to include the dependence of Si on a, which we deduce from Fig. |. 

For the present purposes it suffices to set Si = AD and dSi/da = J2 for all diagrams of 
Fig. |16|, except the last one, for which dSi/da = (J2 — Ji). Thus for the first diagram of 



Fig. |16| we have 

J2 <^TlEz{<y, V) = -l2JlJl{Ab)~\dS/da) . (A7) 

Indicating with 5Vi{2, 2) and 5Vi{2, 3) the total contribution to ¥^{2, 2) and ¥^{2, 3) from 
the diagrams of Fig. |16| one has 

51/4(2, 2) = -{Jl/D^) (aSJI + 24 J2 _ 12 Ji J2) , (A8) 
514(2,3) = -(6 Jl/^^) . (A9) 

However, one also needs to consider terms proportional to where two pairs of excita- 
tions are created and destroyed which do indeed turn out to be important. Consider first a 
set of four spins rii at sites i and the following processes 
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(i) ni,n2 excited, ni,n2 destroyed, n3,n4 excited, n3,n4 destroyed 

(ii) ni,n2 excited, n'^.n^ excited, n'^.n^ destroyed, ni,n2 destroyed 

(iii) ni,n2 excited, n3,n4 excited, ni,n2 destroyed, n3,n4 destroyed 

(iv) n^^n^ excited, n3,n4 destroyed, ni,n2 excited, ni,n2 destroyed 

(v) n3,?T.4 excited, ni,n2 excited, ni,n2 destroyed, n3,n4 destroyed 

(vi) 77,3,714 excited, 77.1,712 excited, 713,714 destroyed, 711,712 destroyed 



We will be interested in the cases shown in Fig. where tii and 712 must be first or second 
neighbors straddling one wall and similarly for 713 and 714 with respect to the other wall. 
Except for the possibility that 712 = 713, all the 7i's are distinct. Because we are working to 
linear order in l/S" (that is ignoring terms higher than quadratic in the boson Hamiltonian) 
the energy denominator depends on position on the lattice but not on the position of the 
other excitations. Hence the energy denominators are simply the sum of the energies of 
the excited spins relative to the ground state energy. We denote them by Eijk... when spins 
i,j,k... are excited. Noting that the matrix elements, say M, are common to all processes 
(i)-(vi) we are now in a position to write down the contribution from these diagrams to 
fourth order in perturbation theory 

^ f / I \ ^ ^ ] (AlO) 

S V -^12-^34 -£'12-£'34 E12E1234E12 £'i2-El234-£'34 E34E1234E34 J 

where ^ = 2 if two excitations are present at the same site (Bose statistics) and ^ = 1 for 
all spins distinct. 

Using E1234 = E12 + Eu 

Vi = |(l-0%^. (All) 

-^12-^34 

Putting ^ = 1 it is immediately apparent that there is no contribution from diagrams for 
which all Tij are different. There are however terms 0{1/D^) when ^ = 2. Diagrams of 
this type which contribute to V3(2, 2) are shown in Fig. ^a. Only terms with rj = 1 give a 
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contribution different from Eq. Therefore when the sum over a is taken the term proportional 
to cr is the lowest order which survives. Including a factor 2 for diagrams symmetric with 
respect to reflection in the center wall of Fig. ^a one obtains 



Vr (2, 2) = ( 2 J| - J,J^ + 2Ji4 ] , (A12) 



where the superscript indicates a contribution of type a in Fig. |T7 



Similarly the contributions of this type to V^{2, 3) is shown in Fig. |l7| c. They give 

^W(2,3) = (A13) 

There is one further contribution to 1^(2,2). Consider the following order of excitation 
of four spins 

(i) ni,n2 excited, n3,n4 excited, n2,'^^4 destroyed, ni,n3 destroyed 

(ii) ni,n2 excited, ns,n4^ excited, ni,n3 destroyed, n2,n4 destroyed 

(iii) 713,77,4 excited, ni,n2 excited, ni,?7,3 destroyed, 712,77.4 destroyed 

(iv) ns,n4 excited, 711,712 excited, 712,714 destroyed, 711,713 destroyed 

The pairs (711,712), (713,714), (711,713), (712,714), must all be first or second neighbors spanning 
a wall. This means that the only contribution of this type is to V4(2, 2) and is shown in Fig. 



17b. Proceeding as before the sum of all orderings gives 



^ E1234 ^^^^^ 

S i?i3i?24-El2-E34 

where we have included a factor 2 for the reverse order of the perturbations. Evaluating this 
for the relevant diagram 

where a factor 2 for the mirror image process has been included. 
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Finally obtaining 1^3(2,2) from Eqs. (P), (|M]), (0), (|XT2|) ,and (|XT5|) and ¥3(2,3) 
from Eqs. (|X9|), ([A 21), (|X5|),and ( [A13| ) we are in a position to calculate the sum of the 



contributions to F{2, 3) from all the diagrams of Fig. O. We obtain 



4Q2 



(A16) 



where we have used J2 = Ji/2 + 0{1/DS). Combining this with Eqns. ( |A^) and (|A.y|) we 

get 



F(2,3) 



392J| 



This is negative showing that the (2) : (3) phase boundary is indeed stable. 



(A17) 



APPENDIX B: EVALUATION OF G^(/, J) 

1. Formulation for G 

The best way to get G(i'")(E), as defined in Eq. (41), is as the inverse of the matrix 
Qn = [EI — SH'^''"'^]^^ , where X is the unit matrix. We define 

E + 2D + 2Ji-2J2 = A , E + 2D + 2Ji - J2 = B E + 2D + Ji - J2 = C . (Bl) 

Then for n = 13 the matrix Q„ is 
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(B2) 



and 



Q(l,n) 



{Qu) 



-1 



(-l)^+^-iV„(^,j)/A, 



(B3) 



where Nn{i,j) is the minor of Q„ and A„ = DetQ„. Since from inspection of Eq. (|5TD 
An is thus a common factor in V2{n), it is not necessary to carry out a complete evaluation 
of this quantity. The N's are fortunately much easier to calculate. The first step, is to 
expand N by minors to get it in terms of a matrix in which end effects, embodied in B and 
C, have been eliminated. Thereby we obtain [writing G for G^^'^\E)] 



n~l 



n-1 



C dn-3 + 2Ct/it/2(i„_4 + (2C J2 + Jl J2)dn-5 ~^ 2J^iJ2(i„_g + 7 (B4a 



CJidn-3 + {JIJ2 + BCJ2)dn-i + (Jl j| + BJlJ^ + CJiJ^)d, 



n—5 



+ {CJl + B4 + JtJt)dn-& + 2Ji4dn-7 + Jldn-S (B4b) 

(-l)"+^A„Gi,„ = Jldn-3 + 2BJiJ2dn-, + {2JlJl + B^Jl)dn-^ + {2BJi4 + 2Ji4)rfn-6 

+ (25 J| + JlJt)dn-7 + 2Ji4dn-S + Jldn-9 , (B4c) 

where dn denotes the determinant of the n x n matrix D„ which has no end effects: 
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n 
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n 


n 


•J2 


T 

•Jl 


A 


T 

•Jl 


'J2 




















J2 


-Jl 


A 


-Jl 























J2 


-Jl 
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J2 


-Jl 



(B5) 



Prom here on we will set E ^ 2D + 2Ji - 2J2 = 2D'. 



2. Evaluation of dn 

Expanding the determinant of D about its upper left corner, we have 

dp = — J\dp—\ — J2Adp—2 — JlJ2dp—3 — 1/2(^^—4 • 

We define 

do^l 



(B6) 



(B7) 



and dp — Oioicp< 0. The solution of the recursion relation is obtained from the result that 



F{x) = Yl ^p^^ = ( 1 + Jix + J2Ax'^ + JiJlx'^ + J^x'^ 



p=0 



(B8) 



which, when expanded in powers of x allows us to get dp. For this purpose we write 



F{x)^Jl[l-xJ2yi 

i=l ^ 

where yj is obtained as the solution to 

yf + {Ji/J2)y! + {A/J2)y^ + {Ji/J2)yi + 1 = 



(B9) 



(BIO) 
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We write the four roots of this equation as 



yl = ^ye'^ y2 = -iye-'\ = \ j = -ly-^ e-^\ = 1/^2 = (Bll) 



where y and 8 are positive, with y of order \j Aj J2 ^ 1 and 8 ~ \ jy. Equation ( [B10| ) can 
be solved as a quadratic equation for [yl + 1) from which yi can then be obtained. Then 
y and 5 can be obtained as accurately as needed. 

We now obtain dn by expanding F[x) in powers of x. To do that write 



>=0 ^ ^r=0 



(B12) 



where 



Sr, = y\^yl ^1/2 • • • + yl 



p+i p+i 

1/1 - Z/2 



y^Qv 



(B13) 



and 



J. 



111 1 

p ~^ p— 1 p 

.1/1 1/1 1/2 1/2, 



— ^ = — = 1/"^Q 
I/1I/2 1/^^ 



(B14) 



where 



= ivi^v^ - e*(f-2)^ . . . + (-l)Pe-*P^ 



> ^NP/2 C0S[(P+ 1)<^] 

cos 5 ' 

._.N(p+i)/2^i*±2M 

^ ^ cos5 



p even ; 
, p odd 



(B15) 



So 



C^n — '-^2 1/ ^ ^ Qn—rQr-y 
r=0 



-2r 



To leading order in J2/A, this gives 



(B16) 



dp = J2i/^(— 1)^^^ cos(p5) , p even 



j2V(_l)(p+i)/2 sin(p5) 



p odd 



(BIT) 
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3. Leading Evaluation of V2{n) 



To leading order in (J2/-D') Eq. (|B4a|) gives 



G'2,n-i = i-ir^'^-^d^-s . (B18) 



Asymptotically, 



A„~ J2"A"~ (^4Z}'-^^i^tA...j ^ (4D')"(1 - a)^ (B19) 
where a is the correction of order {J/ D'Y. In fact, to leading order in {J2/D') we have 

G2,n-i = ^^(4D'J2)("-3)/^(-l)^ sin[(n - 2)6] 

= ^^5!:^ f J^y^' , (B20) 

(1 - ao)" ^ ^ ^ ^ ' ^ ^ 

for n even, where = 5(J2/4D')^ and 

G2,n-1 = f (-1)^""^^/^ -^M) • (B21) 



for n odd. This result agrees with Eq. ([58| ) to leading order in J2/D'. For small n these 
results reduce to those given in Sec. III. 

This evaluation is clearly not precise enough to tell whether V2(n), Eq. (|53|), is 
nonnegative. Obviously, when 2n6/7c is close to an integer [or more precisely, when 
G'2,n-i — 2Gi „_i = 0], it is necessary to retain the first higher-order terms which are nonzero 
there. For this purpose we need to keep all the terms in Eqs. ( |5B| ) and (]^7\). Also in evalu- 
ating the Green's functions we have to keep a sufficient number of correction terms in Eq. 
( P16|) . The algebra required for this analysis is too involved to be worth presenting. Instead, 
we give the final result in Eq. (^8]) and discuss how we have numerically verified it. 



4. Numerical Verification of Eq. (54) 



Here we describe the comparison between the analytic results of Eq. (^) and our numer- 
ical evaluation. To clarify the comparison we will work with rather large values of D' / J2- 
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The numerical evaluation was done as follows. In all cases we set Ji = 2J2. First we ob- 
tained y and 5 exactly by solving Eq. ( |B1CI| ). Then we constructed the dn using Eq. ( |B16D . 
We checked that the dn so obtained did satisfy the recursion relation of Eq. (IBS) to one 



part in 10^°. Next we constructed the quantities A„Gjj using Eqs. (B4). To calculate V2{n) 
according to Eq. we set A„ = J2 A" and used the approximation A ~ Aq, where 

AD' ( 5J2 \ 

This approximation only affects slightly the scale of V2(n) because all the G's appearing in 
Eq. (155D are proportional to this factor. 

To numerically verify Eq. (|58D we chose to study n odd, because this case is the first 
where V2{n) approaches zero. In particular we will explicitly discuss only the case of n = 



63 and D' / J2 = 400 for which n\J J2/ {AD) ^ ti/2. For the parameters we used, V2(n) 
became very small. This is because to check the asymptotic forms it is convenient to assign 
values to the parameters well outside anything one would encounter experimentally. For 
instance, V2{n)/J2 became of order 10^^°'^. Obviously, to interpret the numerical results, 
it is convenient to consider quantities in which the exponential decay [A~" in Eq. is 
removed. Accordingly, we list in Table 1 the values of 

^exact ^ jl-^n^lX--{G2,n-l - 2Gi,„_i)' , (B23) 

for a few representative cases where F^^^^^ is close to zero. The subscript 'exact' indicates 
that we evaluated this quantity to double-precision accuracy. (Recall that our precise eval- 
uation is for A^G and not for G itself.) We can write Fi as 

^exact ^ ^5A2Ao"J2"2"-V^ = -^{X / X.T {SVaX" / J) . (B24) 

Since A ~ Aq and KiA" is of order unity, Ff^^^ will be of order unity as desired. According 



to Eq. whose validity we wish to verify, we have the analytic result for Fi. 



^analytic _ / A \ AD^ ^^^^^^^ ^ ^^3^ ^ AD^ ^^^3^^^ ^ ^^3^ ^^^^^ 

VAo/ J2 J2 
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This analytic result predicts that Fi should become zero when S = So, where uSq + ASq = 7r/2, 
from which we get = 0.0249322909. From the solution to Eq. ( |B10| ) we have that 



= 5 + (5/6)53 , ^326) 

which gives A/J2 = 1607.037 or, since A = AD + 4J2, -D/J2 = 400.759, in very precise 
agreement with the the location of the zero of Va from the numerical evaluation given in 
Table 1. So we have confirmed that that Fi is proportional to co^{n5 + 45^). 

Now we turn to the evaluation of VB{n). We are especially interested in Vb at the point 
when VA{n) = 0. To study this quantity we have listed in Table 1 values of 

TTiexact — j2—2n\2\~nfA^ ^ c;/^2 \ 

^2 = -'2 ^n^O y^^l,n^l,n-l — ^^l,n) 

A numerical evaluation to double precision accuracy at the point where Va = is given in 
Table 1 as 

^exact ^ _i_oo4 X 10"^ . (B28) 

The analog of Eq. ( P25|) is 

^analytic ^ _ (3/2) ( A/ Aq)" ( J2/^) ' • (B29) 

In contrast to the previous check, here we actually need to rely on the approximations for 
the various scale factors. Assuming A ~ Aq, the above equation the above equation gives 
^analytic _ _g_g ^ IQ"^, Compared to the exact result. This comparison may not seem 
impressive. However, it does indicate that we have assigned the correct order in J2/D to 
Vb- One notices that F2 does depend on D. An empirical fit to the numerical data indicates 
that a better approximation for Vb (when ~ 0) is to write 

6D 



B 



sx 

Therefore we replaced Eq. ( P29|) by 



+ 14sin(2n5 + 85=^) (^i' ' 



(B30) 
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^analytic ^ ^J^^' [1^^) ^ ^^^^^^'^ "^^^^''^ ^ ^^^^ ' ^^^^^ 

Note that Eqs. (|B29|) and (P31| ) give the same result when Fi = because then sm{2n6 + 
85^) = 0. However, Eq. (|B29| ) reproduces the variation in F2 when Fi is not exactly zero. 

To see an overall comparison between numerical and analytic results we also tabulate 
the ratios 

^. _ ^exact /^analytic ^ ^g32) 

where we used Eq. (B31) for j?^^'^^^^^'^ ^ n jg striking that although the quantities F^^^'^^ and 
^exact ygj.y Qygj- jxiauy decadcs, that the ratios Ri between numerical and analytic results are 
essentially constant, except very near where the quantities pass through zero and a small 
error in the phase shift can cause large variations in these ratios. From this discussion we 
conclude that our numerical results corroborate the complicated algebra leading to Eq. (p8|). 
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TABLES 

TABLE L Numerical evaluations of Vo. Here F = 



Eqs. (§2|) and (|B27D . Also Ri is defined in Eq. (§3|). 
the other data is for n5 ~ 7r/2. 



-Fi + -^2; where Fi and -F2 are defined in 
The data for n = 95 is for n5 ^ 3tt/2. All 



D/J2 

99.000 

100.000 

100.100 

100.180 

100.190 

100.200 

100.205 

100.206 

100.207 

100.208 

100.209 

100.210 

100.211 

100.212 

100.213 

100.214 

100.215 

100.220 

100.300 

100.500 

101.000 



n 
95 
95 
95 
95 
95 
95 
95 
95 
95 
95 
95 
95 
95 
95 
95 
95 
95 
95 
95 
95 
95 



F 

.5223948307 
.0165890715 
.0047456479 
.0002864392 
.0000419448 
.0001330808 
.0001945491 
.0002047595 
.0002142755 
.0002230972 
.0002312245 
.0002386576 
.0002453963 
.0002514408 
.0002567910 
.0002614470 
.0002654088 
.0002748052 
.0019330782 
.0268205406 
.2091773178 



^cxact 

.5134465900 
.0152594233 
.0041642309 
.0003018071 
.0001317981 
.0000312329 
.0000069855 
.0000042185 
.0000021456 
.0000007668 
.0000000821 
.0000000915 
.0000007949 
.0000021923 
.0000042838 
.0000070692 
.0000105486 
.0000383541 
.0028406168 
.0292069915 
.2152168605 



^exact 

.0089482407 
.0013296482 
.0005814170 
.0000153679 
.0000898533 
.0001643137 
.0002015345 
.0002089779 
.0002164211 
.0002238640 
.0002313067 
.0002387491 
.0002461912 
.0002536331 
.0002610748 
.0002685162 
.0002759573 
.0003131593 
.0009075386 
.0023864509 
.0060395426 



^1 
1.611 
1.603 
1.602 
1.599 
1.598 
1.593 
1.583 
1.578 
1.568 
1.546 
1.440 
1.782 
1.660 
1.636 
1.626 
1.621 
1.617 
1.610 
1.602 
1.600 
1.596 



R2 
.759 
.705 
.637 
-.112 
-2.200 
2.939 
1.933 
1.834 
1.751 
1.680 
1.618 
1.565 
1.517 
1.475 
1.438 
1.404 
1.373 
1.256 
.889 
.811 
.786 



765.000 



87 



.0009103506 



.0008733126 



.0000370379 



1.058 



1.429 
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765.499 


87 


-.0000018986 


.0000002010 


-.0000020996 


1.057 


.988 


765.500 


87 


-.0000020258 


.0000001521 


-.0000021779 


1.057 


.999 


765.501 


87 


-.0000021463 


.0000001100 


-.0000022563 


1.057 


1.009 


765.502 


87 


-.0000022599 


.0000000747 


-.0000023346 


1.056 


1.018 


766.000 


87 


.0007858863 


.0008272135 


-.0000413271 


1.058 


1.367 


508.000 


71 


.0099390222 


.0097440109 


.0001950113 


1.072 


1.155 


509.000 


71 


.0007640012 


.0007157580 


.0000482432 


1.072 


1.165 


509.345 


71 


.0000015730 


.0000037895 


-.0000022165 


1.071 


.947 


509.350 


71 


-.0000004279 


.0000025192 


-.0000029471 


1.071 


.991 


509.378 


71 


-.0000068571 


.0000001812 


-.0000070383 


1.074 


1.080 


509.379 


71 


-.0000069369 


.0000002476 


-.0000071844 


1.074 


1.081 


509.380 


71 


-.0000070063 


.0000003243 


-.0000073305 


1.074 


1.082 


509.381 


71 


-.0000070653 


.0000004113 


-.0000074766 


1.073 


1.084 


509.400 


71 


-.0000062236 


.0000040289 


-.0000102525 


1.072 


1.102 


509.450 


71 


.0000138198 


.0000313758 


-.0000175559 


1.072 


1.123 


509.600 


71 


.0002289421 


.0002683972 


-.0000394551 


1.072 


1.140 


510.000 


71 


.0019384519 


.0020362219 


-.0000977700 


1.072 


1.149 


390.000 


63 


.7804778656 


.7781343381 


.0023435276 


1.083 


1.021 


395.000 


63 


.2225976723 


.2213686283 


.0012290440 


1.082 


1.026 


400.000 


63 


.0039724306 


.0038217672 


.0001506633 


1.081 


1.029 


400.500 


63 


.0004905488 


.0004457817 


.0000447672 


1.081 


1.023 


400.600 


63 


.0001919976 


.0001683678 


.0000236298 


1.081 


1.014 


400.700 


63 


.0000259217 


.0000234154 


.0000025063 


1.081 


.881 


400.710 


63 


.0000165981 


.0000162033 


.0000003948 


1.081 


.493 


400.720 


63 


.0000085986 


.0000103153 


-.0000017167 


1.081 


1.379 


400.730 


63 


.0000019232 


.0000057512 


-.0000038280 


1.080 


1.164 


400.740 


63 


-.0000034282 


.0000025109 


-.0000059391 


1.080 


1.114 
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400.750 
400.757 
400.758 
400.759 
400.760 
400.761 
400.762 
400.770 
400.780 
400.790 
400.800 
400.900 
401.000 
402.000 
403.000 
404.000 
405.000 
410.000 



63 
63 
63 
63 
63 
63 
63 
63 
63 
63 
63 
63 
63 
63 
63 
63 
63 
63 



.0000074557 
.0000094872 
.0000097244 
.0000099484 
.0000101592 
.0000103568 
.0000105411 
.0000115389 
.0000115948 
.0000103270 
.0000077355 
.0000909688 
.0003219780 
.0098963366 
.0326387797 
.0684928015 
.1174021181 
.5558237738 



.0000005945 
.0000000406 
.0000000144 
.0000000015 
.0000000018 
.0000000153 
.0000000421 
.0000007328 
.0000027875 
.0000061658 
.0000108676 
.0001306676 
.0003827584 
.0101671709 
.0331182868 
.0691796094 
.1182948632 
.5577260578 



.0000080501 
.0000095278 
.0000097389 
.0000099499 
.0000101610 
.0000103721 
.0000105832 
.0000122718 
.0000143824 
.0000164928 
.0000186032 
.0000396988 
.0000607804 
.0002708343 
.0004795072 
.0006868078 
.0008927450 
.0019022840 



1.079 
1.072 
1.067 
1.037 
1.123 
1.095 
1.090 
1.083 
1.082 
1.082 
1.082 
1.081 
1.081 
1.081 
1.081 
1.080 
1.080 
1.079 



1.091 
1.082 
1.081 
1.080 
1.079 
1.078 
1.077 
1.070 
1.065 
1.061 
1.057 
1.044 
1.040 
1.036 
1.036 
1.037 
1.038 
1.043 
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APPENDIX C: ANALYSIS OF EQ. (HD WHEN ONE TERM DOMINATES 



Here we analyze Eq. ( j53| ) in the limit when n is so large that only a small range 
of a is important. Superficially, if only a single value of a were important, one could 
make V2{n) negative by adjusting J2/D so that the first square bracket vanished. We now 
show that this reasoning is incorrect. The argument is most easily described when one 
arbitrarily sets 0a(O) = 1 and 1) = 0. Since these quantities depend only weakly on 
a, this simplification is only a matter of convenience. The crucial a-dependence is that 
m G'„(e(7°°'0)) when n is arbitrarily large. For that regime we treat the case when the 
contribution to V2(n) comes from near where the summand is maximal. Therefore we have 
[writing ta for e^~°°'°)] 



= e^o + 47(a - ao)^ = 2D" + 470;^ 



(CI) 



where D" is of order D and x = a — ao- Thus, when the sum over a is replaced by an 
integral, Eq. (|53D (for n even) becomes 

J2 ' 



Voin) 



dx 



sm'[nJj2/{4D + 47^2)] - C{J2/D) 



(C2) 



AD + 47x2^ 

where we used the result of Eq. ( ^Sf ) to write the negative correction term. Also from 
Appendix B we identified AD' as begin + 2D' and we dropped the distinction between D' 
and D". For large n we write 



\adJ 



1 



2C{J2/D)^^(^^^ ' -Reje^v^ y~ dxe-"^""'+''"'- j 



, (C3) 



where e = 7/-D, e = e — {i'^ /2)^ J2I D'^ , and rj = {?>/ A)in'~^'^ ^ J2I D^ . For large n we can drop 
the term in 77, so that 



V2{n) 



(-) 

\ad) 



TT 



line) 



1 - 2C{J2/Df - Rei e*"v^^ 



1 _ i hh. 
^ 2y D 



(C4) 
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Thus 



2V D I 



Thus V2{n) is positive in the hmit of asymptotically large n. It is easy to see from Eq. (|C2|) 



that for large n the variable x can be of order JD/ (727), which can cause a variation in the 



argument of the sine function of order y J2/-D. This estimate immediately explains the final 
result. 
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FIGURES 

h 




1/D V 

FIG. 1. Schematic representation of interface layering transitions in the Heisenberg model with 
strong uniaxial spin anisotropy, as the magnetic field h changes sign. All interface phases k 
appear in the phase diagram. 



+ + 

i = 1 2 3 k-2 k-1 k k+1 k+2 

FIG. 2. The process which gives the lowest order contribution to the energy difference 
AE'j. = — Ek_i between the interface at positions k and k — 1. + (-) denotes the creation 
(destruction) of a spin excitation by V^. An arrow is used to denote a hop mediated by V||. The 
process shown contributes to E^. but not to Ek_i because the i = 1 spin can not be flipped when 
H ^oo. 
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a a 



Ti ri 



-Tl -Tl 



-c n * -t — m » 

FIG. 3. Configurations needed to calculate the interaction energy for two walls at separation 
n (top) and three walls at separation n and m (bottom). When a = +1 (r) = +1) the left-most 
(right-most) wall is positioned as shown. When a = —1 (77 = —1) the left-most (right-most) wall 
does not exist. 



+ + 



FIG. 4. The contribution from second order perturbation theory which renormalizes the wall 
energy. -|- (-) denotes the creation (destruction) of a spin excitation by Vy/. 



t t t 



^ ^ ^ ^ ^ 



t t 
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FIG. 5. Contribution to the two-wall interaction V2(n) for n = 5 in analogy with the unbinding 
problem of Fig. |2|. This process contributes to V2(5) at order JI/Z)^. 



E-2D 



2J, 



WW 
WW 

oW 



^ 

FIG. 6. The energy, E, of an excitation as a function of position near a wall. We give E — 2D 
when the excitation is created at the circled site. Thus when the excitation is on a site next nearest 
neighboring to the wall its energy is ei = 2D + 2Ji — J2(l — rj) ™d when nearest neighboring its 
energy is 62 = 2D + ( Ji — J2) (1 — ^) • When the wall is absent, these formula give the correct energy 
of an excitation when it is not near a wall. Thus dei/dr] = J2 and de2/dr] = J2 — Ji- 
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FIG. 7. Excitations which contribute to the 2- wall interaction V2{n) for (a) n odd and (b) n 
even. + (-) denotes the creation (destruction) of a spin excitation by Vy/. An arrow denotes a hop 
mediated by Vii. 
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FIG. 8. Schematic phase diagram of the Heisenberg version of the ANNNI model with strong 
uniaxial spin anisotropy D showing the effect of quantum fluctuations. 
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FIG. 9. The diagrams needed to calculate F(2n — 1, 2n) to leading order: contributions to (a) 
V-i{2n -l,2n- 1) (b) V3(2n, 2n + 1). 
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FIG. 10. The diagrams needed to calculate F{2n, 2n + 1) to leading order: contributions to (a) 
F3(2n,2n) (b) V3i2n,2n + 1). 
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FIG. 11. Contribution to ' in second-order perturbation theory. Top: real space represen- 
tation showing sites i and j near the wall where the excitations are created and sites k and I where 
they are destroyed. Bottom: Feynmann diagram representation. 




FIG. 12. Contribution to the energy of a two-wall configuration from third-order perturbation 

theory. In the real space representation excitations are created at sites i and j and ultimately 
destroyed at sites I and m all near the left wall. 
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FIG. 13. Contributions to 6V as per Eq. ( pq) in which the excitation to the left of the wall 
propagates. Nonpropagating contributions are similar to those in Fig. |^. Note that only the term 
in which i and / are nearest neighbors actually appears. Each diagram occurs four times: twice at 
one wall by interchanging creation and annihilation sites, and twice for interchanging the roles of 
the two walls. 



+ 




FIG. 14. Leading contribution to V2{n) for a three-dimensional system. Here r± and s± denote 
vectors in the plane of the wall. 
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FIG. 15. Example of a term contributing to 1^(2, 2) in third order perturbation theory. 



1 t 



- r| - Ti 



I i 



ftll 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 
+ 



+ 



53 



FIG. 16. Processes which can not occur when the wall is as shown and which therefore carry 
a factor Sj^-i (when the wall is absent 77 = —1 and these processes are allowed). The first six 
diagrams contribute to 1/3(2,2) and the last one to V3(2,3). In the last diagram the right-hand 
block contains three down spins. 
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FIG. 17. Terms which contribute to (a) '{2,2), (b) Vf'{2,2), (c) ^(2,2). The figures 
indicate which spins are excited. The way in which all possible orderings of the excitations are 
accounted for is described in the text (see equations (A8) and (A9) for diagrams (a) and (c) and 
equation (AlO) for diagram (b)). In cases (a) and (b) the diagrams which are mirror images in the 
center wall must also be accounted for by including a factor of 2. 
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